Two-stage ordering of spins in dipolar spin ice on kagome 
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Spin ice, a peculiar thermal state of a frustrated ferromagnet on the pyrochlore lattice, has a finite 
entropy density and excitations carrying magnetic charge. By combining analytical arguments and 
Monte Carlo simulations, we show that a similar magnet on kagome has two distinct spin-ice phases. 
The lower-temperature spin ice is distinguished by an ordering of magnetic charges. We comment 
on possible realizations of this model. 



Frustrated magnets [J Q attract attention of both the- 
orists and experimentalists as models of strongly inter- 
acting systems with unusual ground states and elemen- 
tary excitations. One of the recent surprises was a re- 
alization that elementary excitations in spin ice [H, Q| 
are quasiparticles carrying magnetic charge [Bj]. Sub- 
sequent work [6( uncovered signatures of the magnetic 
monopoles in magnetization dynamics of spin ice. Spin 
ice is a frustrated ferromagnet discovered in the py- 
rochlore Ho2Ti 2 07, where magnetic Ho 3+ ions form a 
network of corner-sharing tetrahedra Q. The magnetic 
moments fii = ifj&i are forced to point along a (111) 
axis e, by a strong crystal held. The easy-axis anisotropy 
makes the spins Ising-like, so that a microstate of this 
magnet can be described by Ising variables <7j = ±1. The 
Hamiltonian includes exchange interactions of strength J 
for pairs of nearest neighbors (ij) and dipolar interactions 
between all spins 
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where D = (/io/47r)/i 2 /r 3 n is a characteristic strength 
of dipolar coupling, are spin locations, = (r$ — 



r i)/l r i 



and 



is the distance between nearest 



neighbors. In the absence of dipolar interactions, D = 0, 
and for ferromagnetic exchange, J > 0, the system is 
strongly frustrated because it is impossible to minimize 
the energy of every bond (ij). In a ground state, two 
spins point into every tetrahedron and two point out, 
which is reminiscent of proton positions in water ice Q , 
where every oxygen has two protons nearby and two far- 
ther away. This ice rule is satisfied by a macroscopically 
large number of microstates, so that both protons in wa- 
ter ice and magnetic moments in spin ice can remain 
disordered even at low temperatures 

Large magnetic moments (fj, = 10^s in Ho2Ti 2 07) 
make magnetic dipolar interactions between nearest 



neighbors comparable to exchange [l(| • Together with 
the long-distance nature of dipolar interactions, the sub- 
stantial value of D casts doubt on the usefulness of the 
short-range (D = 0) model of spin ice. Yet numeri- 
cal simulations show that, even after the inclusion of 
dipolar interactions, energy differences between states 
obeying the ice rule remain numerically small — so small 
that magnetic order induced by the dipolar interactions 
is expected to occur only at a rather low temperature, 
T w 0.08D 11J ( we set fee = 1). The apparent weakness 
of the dipolar interaction was explained by Castelnovo et 
al. [5(, who introduced a "dumbbell" version of spin ice, 
in which magnetic dipoles are stretched into bar mag- 
nets of length a such that their poles meet at the centers 
of tetrahedra. The energy of the resulting model can 
be represented as a Coulomb-like interaction of magnetic 
charges at the ends of the dumbbells, qi — ±n/a [5(: 
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In this expression, Q a — J2ie a 1i IS tne sum °f magnetic 
charges at the center of tetrahedron a. In a spin-ice state 
of the dumbbell model, every tetrahedron will have two 
north and two south poles with a total magnetic charge 
Q a = 0, minimizing the first term in Eq. |2|). As a re- 
sult, no magnetic field will be generated and the magnetic 
dipolar energy will be strictly zero. A similar, albeit in- 
complete cancellation occurs in the original model, mak- 
ing Eq. @ a very good approximation for the energy (jTJ) . 
The charge of tetrahedron a, expressed in units of fi/d, 
is calculated as 



(3) 



with the plus sign for one sublattice of tetrahedra and 
minus for the other. Residual interactions, responsible 
for the formation of magnetic order, are weak and fall 
off quickly with the distance [||. The resulting energy 
differences between states obeying the ice rule are only 
a small fraction of the dipolar energy scale D 12|, 1 1 31 ] - 
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FIG. 1: Magnetic configurations of the dipolar kagome ice and 
their alternative representations, (a) A spin-ice microstate 
lacking spin order but possessing charge order. The latter is 
manifest in the dumbbell representation (b). (c) One of the 
ground states exhibiting the y/3 x \/3 magnetic order and its 
depiction in terms of dimers (d). The dashed line marks the 
magnetic unit cell. 



As the magnet is cooled down from a high-temperature 
paramagnetic state with completely uncorrelated spins, it 
first gradually enters the spin-ice regime at the crossover 
temperature T ss 2J e g, where J c g = J/3 + 5D/3 is the 
effective interaction for nearest-neighbor Ising spins Uj 
0, [HI , and then undergoes a phase transition into a mag- 
netically ordered state at T w 0.08ZX 

In this Letter we discuss a related problem of dipolar 
spin ice on kagome, a two-dimensional lattice of corner- 
sharing triangles. Each spin is constrained to point along 
the line connecting the two triangles, Fig.QJa). Possible 
ways of realizing such a system are discussed below. A 
short-range version of this model (D = 0, J > 0) was 
studied by Wills et al. 

At a first glance, the dipolar spin ice on kagome is 
very similar to its counterpart on the pyrochlore lattice 
and one might expect a similar sequence of states: as 
the system cools down, it would gradually enter a corre- 
lated but still disordered spin-ice state and then make a 
phase transition into a magnetically ordered state. How- 
ever, a closer look reveals the existence of an intermediate 
thermodynamic phase with disordered spins and ordered 
magnetic charges. 

The existence of the intermediate charge-ordered phase 
depends on the allowed values of magnetic charge ([3]), 



which can be even on a tetrahedron and odd on a tri- 
angle. Minimization of the first term in Eq. (J2]) yields 
Q a = on the pyrochlore lattice and ±1 on kagome. 
A microstate obeying the ice rule Q a = ±1, shown in 
Figs. [Ha) and (b), contains uncompensated charges that 
generate a magnetic field. To a first approximation, 
the system energy is given by the Coulomb terms (|2|). 
Nonzero values of magnetic charges results in substan- 
tial energy differences between different states obeying 
the ice rule. The Coulomb energy is minimized when the 
charges are ordered in such a way that adjacent trian- 
gles carry charges of opposite signs. The charge-ordered 
states of the dipolar kagome ice are closely related to the 
ice states of the p yro chlore spin ice in a magnetic field 
along a (111) axis [151 ], The number of such states grows 
exponentially with the number of spins N. They are ex- 
actly degenerate in the dumbbell model. In the dipolar 
spin ice, the degeneracy is broken by small corrections to 
the Coulomb energy ([2]). 

This hierarchy of energy scales suggests the following 
sequence of thermal transformations in the dipolar spin 
ice on kagome. As the magnet cools down from the high- 
temperature paramagnetic state with entropy per spin 
s = In 2 = 0.693, it gradually enters the spin-ice state 
with s w (1/3) In (9/2) = 0.501 Q. At a temperature of 
the order of D, it will enter a distinct phase with ordered 
magnetic charges, where entropy density is reduced but 
remains nonzero, s = 0.108 per spin; the state is similar 
to that of the pyrochlore spin ice in a magnetic field along 
(111) [15]. At an even lower temperature, the system will 
enter a spin-ordered state with zero entropy density. In 
contrast, the short-range version of spin ice with nearest- 
neighbor interactions exhibits neither charge, nor spin 
order 14 1. 

The spin order emerging on top of magnetic charge 
order is expected to be the of y/3 x y/3 type shown in 
Fig. HJc), the same as in the short-range model with an- 
tiferromagnetic second-neighbor exchange [3| ■ This can 
be understood as follows. In a charge-ordered state ev- 
ery triangle has two majority spins pointing into (or out 
of) the triangle and a minority spin pointing the other 
way. Such states can be represented by dimer cover- 
ings of a honeycomb lattice, Fig. [ljd); the dimers in- 
dicate locations of minority spins. The energy of such 
a state is determined by the interactions between mi- 
nority spins alone. To see that, picture a minority spin 
—\x as a superposition of a majority spin +fi and a mi- 
nority spin of double strength — 2/x; in this represen- 
tation, majority spins form an inert background. We 
thus arrive at a model of dimers with point dipolcs of 
strength 2/i directed along the dimers and towards tri- 
angles with positive charge. The interaction energy of 
two dimers depends on their mutual position. It is mini- 
mized by increasing the number of second neighbors (dis- 
tance between centers s/3r nn ) and reducing the number 
of third neighbors (2r nn ). The dimer configuration shown 
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FIG. 3: Temperature dependence of the specific heat c(T) 
and entropy per spin s(T) of the dipolar spin ice ([TJ with (a) 
ferromagnetic exchange J = 0.5D and (b) antiferromagnetic 
exchange J = — 2.67D. The linear size of the system is L = 
12. The dashed lines show levels of entropy s — 0.693 (Ising 
paramagnet), 0.501 (spin ice), and 0.108 (charge-ordered spin 
ice) per spin. 



FIG. 2: Monte Carlo simulation of the charge-ordering tran- 
sition in the dumbbell model ([2}. (a) and (b) show the tem- 
perature dependence of the (normalized) staggered charge or- 
der parameter Q and Binder's fourth-order cumulant B4, 
I — (Q 4 )/3(Q 2 ) 2 , respectively, for various system sizes L. (c) 
and (d) show the scaling of specific heat c and charge-order 
susceptibility xq> respectively. 

in Fig. [TJd) optimizes both. It is one of three states re- 
lated to each other by lattice translations. The corre- 
sponding spin order is shown in Fig. [He). There are a 
total of 6 magnetic ground states related to each other 
by lattice translations and time reversal. 

To verify this scenario, we performed Monte Carlo sim- 
ulations using both the original model of dipolar spin ice 
((T|) and its dumbbell approximation ©. In both cases, 
the independent variables were Ising spins related to 
magnetic charge by Eq. ([3]). We use the lattice constant 
as a unit of length, so that the nearest neighbor distance 
is r nn = 1/2 and the dumbbell length is a = l/\/3. 

Simulations of the dumbbell model were performed on 
lattices with periodic boundary conditions and lengths 
up to L = 30, which corresponds to N = 3L 2 = 2700 
sites. Temperature was measured in units of D = 
(/Lto/47r)/Lt 2 /rn n . The charging energy — the first term in 
Eq. ©—had the coupling strength 1/C = 8.264(^ /4tt). 
The long-range Coulomb interaction in Eq. @ was trun- 
cated beyond the distance of 10r nn . The temperature 
dependence of the staggered charge order parameter Q is 
shown in Fig. [2ja) for various system sizes. The existence 
of a continuous charge-ordering transition is evidenced by 
Fig. d^b), where Binder's fourth-order cumulants cross 
at T c w 0.267 D for different system sizes. 

On the basis of the symmetry of the order parame- 
ter, one might expect the charge-ordering transition to 
be in the universality class of the two-dimensional Ising 
model. However, the long-range nature of the dipolar 



interactions may change that. A finite-size analysis re- 
veals a universal scaling behavior of the specific heat c 
and charge susceptibility \Q = ((Q 2 ) — (Q) 2 )/NT, shown 
in Figs. [2jc) and (d). In obtaining the scaled curves, 
we have subtracted a size-dependent background con- 
tribution from the specific- heat. The critical exponents 
extracted from the analysis, a = 0.09, v = 1.13, and 
7 = 1.93, are close, but not equal to those of the two- 
dimensional Ising model. Given the crude handling of 
the long-range interactions, the critical exponents may 
not be accurate. Further work is required in this direc- 
tion. 

Monte Carlo simulations were also performed for the 
original dipolar model fT]) on lattices with similar sizes. 
The specific heat c(T) and entropy per spin s(T) are 
shown in Fig. [3] for ferromagnetic exchange J = 0.5-D, 
and antiferromagnetic J = — 2.67-D. In both cases, a 
broad peak in c(T) signals a crossover from the param- 
agnetic regime to the spin-ice state. The latter is seen 
as a washed-out plateau in s(T) near the characteris- 
tic spin-ice value s ~ 0.501. The crossover tempera- 
ture T w 2J c g, where J e g = J/2 + 7D/4 on kagome. 
Two sharp peaks in c(T) at lower temperatures mark the 
charge and spin-ordering phase transitions. In the anti- 
ferromagnetic case, J = — 2.67D, the effective nearest- 
neighbor Ising coupling becomes small, J c g — 0.415Z?, 
so that the system is near the spin ice-antiferromagnet 
boundary, given approximately by the condition J e ff = 
This makes the spin- ice plateau indistinct, but the 
two phase transitions are clearly present. Detailed tem- 
perature scans were performed in the vicinity of the mag- 
netic ordering. It is known from previous studies of mag- 
netic ordering in the pyrochlore spin ice [H, [3] that 
simple local Monte Carlo algorithms fail to find mag- 
netic order because a single spin flip violates the spin-ice 
rule, leading to a large increase in the energy and a likely 
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FIG. 4: Monte Carlo simulations of magnetic phase tran- 
sitions, (a) and (b) show the temperature dependence of 
the (normalized) magnetic order parameter M and Binder's 
fourth-order cumulant B 4 = 1 - (| Af| 4 )/3(|M | 2 } 2 , respec- 
tively, for various system sizes L. (c) and (d) show the data 
collapsing of scaled specific heat c and magnetic susceptibility 
X m = «|A/| 2 ) - \(M)\ 2 )/NT, respectively. 



rejection of the new state. A similar problem arises in 
the charge-ordered state of kagome spin ice. To circum- 
vent it, we added nonlocal updates — flipping spins along 
closed loops — that keep the magnetic charges unchanged. 

The results for dipolar spin ice on kagome with antifer- 
romagnetic exchange J — —2.67D are shown in Fig. 2] 
The magnetic order parameter M = (1/iV) J2< a i e jkri 
with k = (47r/3, 0) is shown as a function of temperature 
in Fig. [UJa) for various system sizes. Binder's cumulants 
cross roughly at T m w 0.089 D [Fig. EJb)], indicating a 
continuous transition. 

The symmetry-breaking pattern of the magnetic or- 
der described above suggests that the magnetic transi- 
tion is in the universality class of two-dimensional three- 
state Potts model. Finite-size scaling of specific heat 
and magnetic susceptibility yielded the critical exponents 
a = 0.31, v = 1.09, and 7 = 2.19 for the magnetic 
phase transition. While the exponent a is close to the 
expected exact value 1/3, the other two differ signifi- 
cantly from their counterparts in the two-dimensional 
three-state Potts model (y = 5/6 and 7 = 13/9) 
The discrepancy may be due to the long-range nature of 
dipolar interactions. 

We presented a plausible scenario of a two-stage mag- 
netic ordering in the dipolar spin ice on kagome. The 
system has two spin-ice regimes with strongly fluctu- 
ating spins, distinguished by the ordering of magnetic 
charges in the low-temperature ice phase. A separate 



spin-ordered phase exists at the lowest temperatures. 

Spin ice on kagome can be created by substituting 1 /4 
of magnetic sites of a pyrochlore spin ice with nonmag- 
netic atoms, as was done in the Heisenberg antiferromag- 
net herbertsmithite ZnCu3(OH)gCl2 [17]. Kagome spin 
ice already exists as an artificial magnetic lattice [Lsl . 
Although the energy scale of dipolar interactions in ar- 
tificial ice greatly exceeds the room temperature [20j, it 
may be possible to introduce thermal motion of spins by 
agitating them in a manner of granular systems [211 ]. 
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